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Physical Interpretations of Negative Imaginary Systems Theory 


Ian R. Petersen 


Abstract — This paper presents some physical interpretations 
of recent stability results on the feedback interconnection 
of negative Imaginary systems. These interpretations involve 
spring mass damper systems coupled together by springs or 
RLC electrical networks coupled together via inductors or 
capacitors. 

I. Introduction 

In recent years the theory of negative imaginary systems 
has emerged as a useful complement to positive real theory 
and passivity theory; e.g., see [l]-[3]. It is well known that 
linear mechanical systems with force inputs and collocated 
velocity outputs lead to positive real transfer functions; e.g., 
see [4]. Similarly, linear mechanical systems with force 
inputs and collocated position outputs lead to negative imag¬ 
inary transfer functions; e.g., see [2]. In this paper we extend 
these ideas to consider mechanical interpretations of the main 
negative imaginary stability results and compare them with 
corresponding interpretations of positive real stability results. 
The negative imaginary stability results involve the positive 
feedback interconnection between two negative imaginary 
system, one of which is regarded as the controller and one 
of which is regarded as the plant; e.g., see [5], [6]. In 
our mechanical interpretation of these results, the plant is 
a spring mass damper system and the controller is another 
spring mass damper system which is coupled to it via a 
spring. This is similar to the behavioural systems theory 
notion of control by interconnection; e.g., see [7]. This paper 
also presents similar RLC electrical circuit interpretations of 
the negative imaginary stability results in which the coupling 
between the plant RLC circuit and the controller RLC circuit 
is via capacitive or inductive coupling. 

The mechanical interpretations presented in this paper may 
help to better understand the negative imaginary stability 
results and motivate extensions to these results. Also, the 
mechanical interpretation of the controller may be useful in 
controller design and tuning. 

IT Negative Imaginary and Positive Real 
Systems Theory 

Negative imaginary and positive real systems. Stability re¬ 
sults about the positive feedback interconnection of negative 
imaginary systems take their simplest form in the case in 
which none of the systems have poles at the origin. Hence, 
we will consider this case first along with a stardard positive 
real stability result. In the sequel, we will look at the more 
general case in which poles at the origin are allowed. 

Definition 1 (See [5]): A square real-rational proper 
transfer function matrix G(s) is termed negative imaginary 
(NI) if 


1) G(s) has no poles at the origin and in 3fi[s] > 0; 

2) j[R{j^)—R*{3^)] > 0 for all w G (0,cx)) except values 
of Lo where jw is a pole of G(s); 

3) If jwo, Wo G (0, oo), is a pole of G(s), it is at most a 
simple pole, and the residue matrix Kq = lims^j^j, (s — 
jwo)jG(s) is positive semidefinite Hermitian. 

Definition 2 (See [1 ]): A square real-rational proper 
transfer function matrix G(s) is termed strictly negative 
imaginary (SNI) if 

1) G(s) has no poles in 5ft[s] > 0; 

2) j[G(ja;) - G*(jw)] > 0 for w G (0, cx)). 

In addition to looking at physical interpretations of nega¬ 
tive imaginary stability results, we will also compare these 
interpretations with physical interpretations of positive real 
stability results. Hence, we introduce the corresponding 
notions of positive real and weakly strictly positive real 
transfer functions. 

Definition 3 (See [8]): A square real-rational transfer 
function matrix G(s) is termed positive real (PR) if 

1) No element of G(s) has a pole in 5R[s] > 0; 

2) G(s) -f G*(s) > 0 for sft[s] > 0. 

Definition 4 (See [4]): A non-zero square real-rational 
transfer function matrix G(s) is weakly strictly positive real 
(WSPR) if 

1) G(s) has no poles in 5i[s] > 0; 

2) G{juj) + G*{juj) > 0 for w G (—oo, oo). 

The following lemma makes clearer the relationship be¬ 
tween negative imaginary and positive real transfer functions. 

Lemma 1 (Theorem 2.7.2 of [8]): Let G(s) be a square 
real-rational transfer function matrix. Then G(s) is positive 
real if and only if 

1) No element of G(s) has a pole in 3?[s] > 0; 

2) G(ja;) -I- G*(ja;) > 0 for all real w except values of to 
where juj is a pole of G(s); 

3) If jwo is a pole of any element of G(s), it is at 

most a simple pole, and the residue matrix, Kq = 
lims_>.j^(, (s — jwo)G(s) in case wq is finite, and AToo — 
lim(^_>oo infinite, is positive semidef¬ 

inite Hermitian. 

Now we recall a stability result for the positive feedback 
interconnection of two negative imaginary systems, denoted 
by [G(s),G(s)], as shown in Figure [T] 

Theorem 1 (See [5]): Given a negative imaginary transfer 
function matrix G(s) and a strictly negative imaginary trans¬ 
fer function matrix G(s) that also satisfy G(oo)G(oo) = 0 
and G(oo) > 0. Then the positive feedback interconnection 
[G(s),G(s)] is internally stable if and only if 

Amax(G(0)G(0)) < 1. 


( 1 ) 




Fig. 1: Positive feedback interconnection 


Also, we recall a stability result for the negative feedback 
interconnection of two positive real systems, denoted by 
[C?(s), —G(s)], as shown in Figure |2] 



Fig. 2; Negative feedback interconnection 


Theorem 2 (See for example Lemma 3.37 in [4]): Given 
a positive real transfer function matrix G{s) and a weakly 
strictly positive real transfer function matrix G{s) such 
that their negative feedback interconnection is well posed. 
Then the negative feedback interconnection [G(s),—G(s)] 
is internally stable. 

Negative imaginary systems with poles at the origin. 

We now recall some recent stability results for the feed¬ 
back interconnection between negative imaginary systems 
containing poles at the origin. This situation may arise in 
control systems in which the plant has free body motion 
or the controller contains integral action. The complexity of 
these results depends on the special cases being considered. 
For the purposes of this paper, we will present two of the 
simplest results. 

We begin with a generalized definition of the negative 
imaginary property. 

Definition 5 (See [6]): A square transfer function matrix 
G(s) is NI if the following conditions are satisfied; 

1) G(s) has no pole in 3fJ[s] > 0. 

2) For all w > 0 such that juj is not a pole of G(s), 
3 (G(jo.) - G(jcc)*) > 0. 


3) If s = jwo with Wo > 0 is a pole of G(s), then it is 
a simple pole and the residue matrix K = lim (s — 

s— >ju)o 

jwo)jG(s) is positive semidefinite. 

4) If s = 0 is a pole of G(s), then lim s^G{s) = 0 for 

S -^0 

all k > 3 and lim s^G(s) is positive semidefinite. 

s —^■O 

Then, we define the following matrices which will be used 
in the stability conditions to be presented; 

Ga = lim s^Gis), Gi = lim s f G(s) - ^ ) . (2) 

s ->0 s ->0 V S‘‘ J 

Roughly speaking, transfer function matrices G(s) with 
only single poles at the origin have G 2 = 0 and transfer 
function matrices with only double poles at the origin have 
Gi =0. 

Theorem 3 (See [6]): Suppose that the transfer function 
matrix G(s) is SNI and the strictly proper transfer function 
matrix G(s) is NI with G 2 = 0 and Gi invertible. Then, the 
closed-loop positive-feedback interconnection [G(s),G(s)] 
is internally stable if and only if 

G(0) < 0. (3) 

Theorem 4 (See [6]): Suppose that the transfer function 
matrix G(s) is SNI and the strictly proper transfer function 
matrix G(s) is NI with Gi = 0 and G 2 > 0. Then, 
the closed-loop positive-feedback interconnection between 
[G(s),G(s)] is internally stable if and only if 

G(0) < 0. (4) 

III. Mass Spring Damper System Interpretations 

We first consider a physical interpretation of Theorem |2] 
for the case in which the controller transfer function G(s) = 
d > 0 is a constant and the plant transfer function G(s) 
is a scalar transfer function corresponding to a spring mass 
damper system in which the plant input is the force applied 
to a given mass and the plant output is the corresponding 
velocity of that mass. The equations of motion for such a 
plant can be written in the form 

Mx + Dx + Kx = Lu, y = 

where M > 0, D > 0, K > 0, x € corresponds to the 
vector of mass positions, u G M corresponds to the force 
input and j/ G R corresponds to the velocity output; e.g., see 
[4], [9]. The corresponding plant transfer function is then 
G(s) = (s^M + sD + iT) ^ L. It is straightforward to 
verify that G(s) is PR and G(s) is WSPR. In this case, the 
negative feedback interconnection between G(s) and G(s) 
can be interpreted as shown in Figure [3 Thus in this case, the 
positive real stability theorem. Theorem |2] can be interpreted 
as saying that the addition of a damper to a spring mass 
system will lead to the overall system being internally stable. 
We can think of the damper as a static velocity feedback 
controller; e.g., see [10]. This interpretation can be extended 
to allow for a more general dynamic controller with a WSPR 
controller transfer function G(s) as illustrated in Figure 0] 
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Fig. 3: Spring mass damper interpretation of the PR stability 
theorem 



Fig. 4: Spring Mass Damper illustration of the PR stability 
theorem with dynamic plant and controller 


From this figure, we can write down the equations of 
motion as follows; 

miXi = F — kixi — dixi, 

1712X2 = -k2X2 - d2X2 - d{x2 - Xi). 

The input to the plant is the force F and the output of 
the plant is the velocity xi. Taking Laplace transforms, we 

obtain the plant transfer function G(s) = -—n—■ 

If we assume mi > 0, di >0 and ki > 0, it is 
straightforward to verify that G{s) is PR. For a negative 
feedback interconnection, the input of the controller is the 
output of the plant xi and the output of the controller is 
minus the input of the plant —F = —d{x 2 — ii). Taking 
Laplace transforms, we obtain the controller transfer function 
G{s) = which is WSPR for m^ > 0, 

d > 0, ^2 > 0 and ^2 > 0. Thus, the PR stability result in 
this case can be interpreted as simply ensuring the stability of 
the complete spring mass system shown in Figured The fact 
that the PR stability result can be used to ensure stability is 
directly related to the fact that the plant subsystem is coupled 
to the controller subsystem via a damper for which the force 
is proportional to velocity. Note that in this example, the 
condition d 2 > 0 is required for the controller transfer 
function G{s) to be WSPR. If d 2 = 0 then G{s) will have 
a pair of complex zeros on the imaginary axis. Also, G{s) 
is never WSPR since it always has a zero at the origin. 

We will now see that the NI stability result. Theorem [T] 
corresponds to the case in which the plant subsystem is 
coupled to the controller subsystem via a spring in which 
the force is proportional to displacement. Indeed, consider 
a physical interpretation of the Theorem [T] for the case in 
which the controller transfer function G(s) = —k < 0 is 


a constant and the plant transfer function G(s) is a scalar 
transfer function corresponding to a spring mass damper 
system in which the plant input is the force applied to a 
given mass and the output is the corresponding displacement 
of that mass. The equations of motion for such a system can 
be written in the form 

Mx + Dx + Kx = Lu, y = F’" X 

where M > 0, D > 0, AT > 0, x € R” corresponds to the 
vector of mass positions, m € R corresponds to the force 
input and y G R corresponds to the displacement output; 
e.g., see [4], [9]. The corresponding transfer function is then 
G(s) = L’^ (^s^M + sD + at) ^ L. It is straightforward to 
verify that G(s) is SNI and G(s) is NI. In this case, the 
positive feedback interconnection between G(s) and G(s) 
can be interpreted as shown in Figure |5] Note that in this 
case, we have reversed the strict and non strict transfer 
functions compared to the PR case. Also, in this case, we 
do not actually need to assume that fc > 0 in order to ensure 
stability. Indeed according to Theorem [T] we only require the 
DC gain condition ([T]i which is G(0)G(0) = —A^ATLfc < 1. 
This will be automatically satisfied if fc > 0. However, it will 
also be satisfied if fc < 0 and —k< That is, the NI 

stability theorem can allow for the situation in which the 
controller corresponds to a negative spring provided it is not 
too large. 

From the above example, we can see that the application 
of the NI stability result Theorem [Harises when the plant is 
coupled to the controller via a spring. This interpretation can 
be extended to allow for a more general dynamic controller 
with an NI transfer function G(s) as illustrated in Figure [h] 
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Fig. 5; Spring mass damper illustration of the NI stability 
theorem 



Fig. 6: Spring mass damper illustration of the NI stability 
theorem with dynamic plant and controller 


From this figure, we can write down the equations of 
motion as follows; 

miaii = F — k^xi — di±i 

7712 X 2 = -k2X2 - d2±2 - k{x2 - Xi). 

The input of the plant is the force F and the output of the 
plant is the displacement xi. Taking Laplace transforms, we 
obtain the plant transfer function G{s) = mis ' ^+dis+ki • 
we assume mi > 0, di > 0 and fci > 0 it is straightforward 
to verify that G{s) is SNI. For a positive feedback inter¬ 
connection, the input of the controller is the output of the 
plant xi and the output of the controller is the input of the 
plant F = k{x 2 — xi). Taking Laplace transforms, we obtain 
the controller transfer function G{s) = -k 
which is NI for m 2 > 0, (i 2 > 0 and k + k 2 >0. Thus, the 
NI stability result in this case can be interpreted as simply 
ensuring the stability of the complete spring mass system 
shown in Figure |6] The fact that the NI stability result can be 
used to ensure stability is directly related to the fact that the 
plant subsystem is coupled to the controller subsystem via a 
spring for which the force is proportional to displacement. 

Note that in this case, since the controller is only required 
to be NI, we can allow d 2 = 0 and also, we can allow k 
or /c 2 to be negative provided that k + k 2 >0 and the DC 
gain condition ([T]i is satisfied. For this example, the DC gain 
condition ([T]i is G(0)G(0) = < 1. If both k and 

^2 are positive, this condition will be automatically satisfied. 
However, if one of k and ^2 is negative, this condition will 
require an additional restriction of In this 

case, this condition will in fact imply that the condition k + 
A :2 > 0 is satisfied. 


A distinction between the use of the PR stability theorem 
and the NI stability theorem is that the PR stability theorem 
requires that all components of the controller to be passive 
whereas the NI stability theorem allows for non-passive 
spring components in the controller. However, since the 
application of the NI stability theorem requires that the plant 
be SNI this means that in this example, the plant damping di 
was required to be positive whereas in the the application of 
the PR stability theorem, the plant damping could be zero. 

We now consider a mass spring damper interpretation 
of Theorems [3 and |4] corresponding to a plant transfer 
function with poles at the origin. In our mass spring damper 
interpretation, this will correspond to a plant with free body 
motion such as shown in Figure [T] 

From this figure, we can write down the equations of 
motion as follows: 


miXi = F — diXi, 

7712 X 2 = —fc2X2 — C?2X2 — fc(x 2 — Xl). 


The input of the plant is the force F and the output of the 
plant is the displacement xi. Taking Laplace transforms, we 
obtain the plant transfer function 


G{s) 


1 

-I- 


(5) 


If we assume ttii > 0 and di > 0 it is straightforward 
to verify that G(s) is NI according to Definition |5] For a 
positive feedback interconnection, the input of the controller 
is the output of the plant xi and the output of the controller 
is the input of the plant F = k{x 2 — xi). Taking Laplace 
transforms, we obtain the controller transfer function G(s) = 





































transfer functions. 



Fig. 7: Spring mass damper illustration of the NI stability 
theorem with free body plant and a dynamic controller 


IV. Electrical Circuit Interpretations 

We now consider electrical circuit interpretations of the NI 
stability result Theorem [T] We first consider an RLC circuit 
interpretation of Theorem [T] as shown in Figure | 8 ] where 
the plant input u corresponds to the charge g = J idt on 
the controller coupling capacitor C and the plant output y 
corresponds to the voltage V on the plant circuit. To find the 



Fig. 8 ; RFC circuit illustration of the NI stability theorem 


-k rTXtnlll. '"hich is SNI for > 0, > 0, A: ^ 0, 

k 2 0 and fc + A :2 > 0 . 

If we assume di > 0, then the conditions of Theorem [3 
will be satisfied provided the DC gain condition Q holds: 


G(0) 


kk2 
k + k2 


< 0 . 


( 6 ) 


This condition will be satisfied if fc > 0 and ^2 > 0 and in 
fact it can only be satisfied if fc > 0 and /c 2 > 0. Thus, in 
this case, we cannot allow the use of active components in 
the controller. 

Note that this example can also illustrate the use of 
Theorem |3 when the controller includes integral action. In 
this case, we reverse the role of the plant and the controller 
in Figure Q In this case, the plant is a spring mass damper 
system which is coupled via a spring such that the input to 
the plant is the displacement of the string and the output of 
the plant is the corresponding force provided by the spring. 
Also, the controller with integral action is given as in Q. 

In the case that di = 0 then the conditions of Theorem 
|4]will be satisfied provided the same DC gain condition (| 6 ]l 
is satisfied. This corresponds to the case in which the plant 
transfer function has a double pole at the origin. 

Note that in this section, we have restricted attention to 
spring mass damper systems operating in one dimension and 
such that all transfer functions are SISO. However, it would 
be straightforward to construct similar examples involving 
spring mass damper systems operating in two or three 
dimensions. Also, it would be straightforward to construct 
examples involving multiple spring couplings between the 
plant and controller which would correspond to MIMO 


plant transfer function, we first write down the impedance of 
the plant circuit: Zi(s) = ^—— 7 ^. Then 


sq 


Hence, the plant transfer function is G{s) = — 2 ^^ 1 1 

S ^1+ + Li 

which is SNI. 

Also, to find the controller transfer function, we write 
down the impedance of the controller circuit: ^ 2 ( 3 ) = 
1 , 1 , ■ Then, 

V 2 = iZ2 = sqZ2 = 


sq 


Rl + + 'SG2 


Hence, q = C{V - V 2 ) = CV- 
CV 


Cs^q 


q = 


1 + 




= c 




T- and therefore 


s2(C' + C 2 ) + ^ + —^ 


-V. 


That is, G(s) = G r - which is NI provided 

C > 0 , G + G 2 > 0 , i ?2 > 0 and L 2 > 0 . 

In this case, the DC gain condition ([Til of Theorem [T] is 
G(0)G(0) = 0 X C < 1 which is automatically satisfied. 
Thus, in this case. Theorem [T] can be interpreted as saying 
that the coupling of the two RFC circuits via a capacitor 
will lead to the overall system remaining internally stable. 
In addition, G(s) will retain the NI property even if G 2 is 
negative provided that C + G 2 >0. Thus, in this case, we 
can also allow for controllers with active elements. 



















































An alternative RLC circuit interpretation of Theorem [T] 
involves inductive coupling rather than capacitive coupling. 
To illustrate this, consider the RLC circuit shown in Figure 
|9l In this circuit, the plant input u corresponds to the flux 
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Fig. 9: RLC circuit illustration of the NI stability theorem 


(j) = J Vdt on the controller coupling inductor and the plant 
output y corresponds to the current ii of the plant circuit. 
To And the plant transfer function, we Arst write down the 
impedance of the plant circuit: Zi(s) = i?i + sLi + 

Then V = —iiZi implies 


V 


s4> 


Zi Ri + sLi + 

Hence, the plant transfer function is G(s) = — ^i 
which is SNI provided Li > 0, i?i > 0 and Ci > 0. 

Also, to And the controller transfer function, we write 
down the impedance of the controller circuit: ^ 2 ( 3 ) = 
i?2 + SL 2 


Then, 

SC 2 

V 




i?2 


*^2 + ^ 


Also, i = Hence, 

’ sL /• ’ 


ii = i + i2= -r + 


Sfj) 


R 2 + SL 2 + ^ 


~{L + L 2 ) + si?2 ^ 


and therefore 6 = 


L{^s^ L2 + SR2 + 

L{s^ L‘2 + sR2 + -^) 


L(s2(L+L2)+si?2 + ■^) 


ii That is, G{s) = 


provided L>0,L + L 2 > 

0, i ?2 > 0 and C 2 > 0. 

In this case, the DC gain condition ([T]i of Theorem [T] is 
G(0)G(0) = 0 X L < 1 

which is automatically satisfled. Thus, in this case. Theorem 
[T]can be interpreted as saying that the coupling of the two 
RLC circuits via an inductor will lead to the overall system 
remaining internally stable. 

V. Conclusions 

In this paper we have presented a number of physical 
interpretations of stability results for feedback interconnec¬ 
tions of negative imaginary systems and compared these with 
stability results for positive real systems. The physical inter¬ 
pretations have involved both spring mass damper systems 
and RLC electrical networks. The results may be useful in 
gaining a better understanding of these stability results and 
motivating new extensions to the existing results. Also, the 
interpretations of a controller as a physical system may be 
useful in robust controller design and tuning even in the case 
when the controller is actually implemented using actuators, 
sensors and computers; e.g, see [11]. 
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